*«»  MOST  BroJect 


NAVAL  UNDERSEA  WARFARE  CENTER 


£ V PUPPED  REFERENCE  CROSS-CORRELATION 
" 1 COMPUTER  PROGRAM . 


W.  J./Wilsterman,  Jr, 


SUBPROJECT  NO|  SF101Q316 
TASK  NO.  11197 


DISTRIBUTION  Si.. 

Approved  lor  putTl-j  r 
Distribution  Unia.i:. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 

A.  Basic  Theory  of  Cross-Correlation 

B.  Clipped  Reference  Cross-Correlation 

C.  Programming  Options 

II.  PROGRAM  USAGE 

A.  Card  Deck  Preparation 

B.  Parameter  Cards  Preparation 

C.  END  Cards 

III.  PROGRAM  OUTLINE:  GENERAL  DESCRIPTION 

A.  CARDINPUT 

B.  TRANTIME 

C.  READPACK 

D.  RDUNPACK 

E.  COMPCORR 

F.  PRINTOUT 

IV.  CARDINPUT 

A.  Translation  of  First  Card  Parameters 

B.  Translation  of  Second  Card  Parameters 

C.  Recognition  of  End  Card;  Exit  Condition 

V.  TRANTIME 

A.  Specific  Explanation 

B.  Compute  REFNUM  and  ECHNUM 


PAGE 

1 

1 

4 

4 

9 

9 

10 

12 

12 

12 

12 

12 

12 

13 

13 

13 

13 

14 

15 
15 
15 
15 


DISTRIBUTION  STM^MUNT  h 

Approved  for  public  roloase; 
Distribution  Unlimited 


t 


TABLE  OF  CONTENTS  (Continued) 


PAGE 


VI.  READPACK  15 

A.  Specific  Information  15 

B.  POSTIME  16 

C.  Mag  Tape  Input  Data  Format  1 6 

D.  Core  Data  Format  18 

E.  SHOVEMIN  19 

VII.  RDUNPACK  19 

A.  PASSFORD  20 

B.  AREA  Table  Design  20 

C.  Output  Function  Words  21 

VIII.  COMFCORR  22 

A.  Continuous  Mode  22 

B.  Segmented  Mode  23 

C.  GTFLAG  Modes  23 

D.  FKFLAG  Modes  24 

E.  UPDATE  S/R  24 

F.  GENCORR  25 

IX.  PRINTOUT  30 

A.  Correlation  Value  Print-Out  30 

B.  Cross-Correlation  Graph  Print-Out  30 

C.  UNPACKTIME  31 

D.  SETEMUP  31 

E.  GRAPH  31 

F.  PRTDATE  32 

G.  SETTIME  32 


ii 


1 


1 « 

I 

I M 

I " 

| 1 

I r 

I 

* < 


TABIiE  OF  CONTENTS  (Continued) 


TIMING 


A.  Data  Input 


B.  Cross-Correlation  Time  Estimation 


C.  CC  Value  Print-Out  Time 


D.  Cross-Correlation  Graph  Print-Out  Time 


P.EFERENCES 


APPENDIX  A 


APPENDIX  B 


Clipped  Reference  Cross-Correlation 
Clipped  Reference  Cross-Correlation 
Clipped  Reference  Cross-Correlation 


Clipped  Reference  Cross-Correlation 


APPENDIX  C 


CLIPCORR  S/R 
READPACK  S/R 


RDUNl’ACK 


COMl'CORR 


PRINTOUT 


PREFACE 


This  memorandum  describes  a computer  program  which  was  developed 
by  the  author  to  analyze  sonar  data  using  cross-correlation  techniques. 
This  program  was  written  for  Code  d603,  Naval  Undersea  Warfare  Center., 

San  Diego  Division,  Subproject  SF  101  03  16,  Task  11197.  This 
memorandum  has  been  prepared  because  it  is  believed  that  the  information 
may  be  useful  in  this  form  to  others  at  the  Naval  Undersea  Warfare  Center 
(NUWC)  and  to  a few  persons  outside  NUWC.  This  memorandum  should  not 
be  construed  as  a report  as  its  only  function  is  to  present  for  the 
information  of  others  an  explanation  of  the  above  mentioned  computer 
program. 


1 


iv 


ABSTRACT 


IRM«m 


n)  The  CLIPREF  System  is  a program  which  performs  clipped 
reference  cross-correlation  between  two  lengths  of  formatted  digitized 
sonar  acoustic  data.  The  program  was  written  by  the  author  for 
Code  D603  in  conjunction  with  the  analysis  of  sea  test  data.  The 
CLIPREF  System  has  a number  of  programming  options  including  selection 
of  data,  segmented  or  continuous  cross-correlation  and  optional  outputs. 
It  also  makes  use  of  certain  time-saving  computation  techniques. 


I.  INTRODUCTION 


A.  Baste  Theory  of  Cross-Correlation 

Definition  of  cross-correlation  function: 

Let  f = f(x)  and  g = g(x)  be  two  integral  real  functions 
over  the  domain  as  xs  b.  Let  t be  a small  variation  of  x.  The 
normalized  cross-correlation  function  a = o(t)  is  defined  as: 
b-r 

f f(x)g(x+r)dx 

o(t)  ■=  a (l) 

b-r 

J |f(x)g(x+T)dx 

a 

From  the  definition  of  the  cross-correlation  function,  it  is  clear 
that  the  range  of  o(t)  lies  between  -1  and  +1.  The  function  is  useful 
in  statistics  because  it  is  a measure  as  to  how  close  two  functions 
resemble  each  other  over  a given  interval.  A high  correlation  value 
of  .8  or  greater  indicates  that  f and  g are  very  similar  over  the 
interval  a%  xsb.  A value  of  +1  indicates  full  positive  correlation; 
both  functions  always  have  the  same  sign.  A high  negative  correlation 
value  of  -.8  or  less  indicates  that  f and  g are  similar  but  opposite 
in  sign.  A correlation  value  of  -1  indicates  that  f and  g are  exactly 
opposite  in  sign  throughout  the  interval  as  x^b.  Low  correlation 
values  from  -0.2  to  +0.2  indicate  that  f and  g bear  little  relation 
to  each  other.  Intermediate  correlation  values  signify  a continuum 
of  evaluations  from  positive  correlation  to  no  correlation  to  negative 
correlation. 

If  f and  g are  distinct  functions  over  as  xs  b,  then  o(t)  is 
called  the  cross-correlation  function  of  f and  g over  the  interval 
as  xs  b.  If  f and  g are  identical  functions,  then  o(t)  becomes 
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the  auto-correlation  function  of  f over  the  interval  as  xs  b. 


Normalized  auto-correlation  is  defined  as: 


b-T 

f f(x)f(x+r)dx 
ff( r)  ■»  a 

J*  |f(x)f(x+T ) |dX 

a 


A digital  computer  cannot  handle  f and  g as  continuous  arbitrary 
real-valued  functions.  It  is  not  possible  to  define  such  a function  in 
advance,  i.e.,  a priori,  in  order  to  program  the  computer  to  perform 
the  integration  and  to  compute  the  function  o(t).  So  in  practice, 
the  function  o(t)  is  estimated.  This  is  done  by  estimating  first 
the  functions  f and  g. 

Consider  the  arbitrary  continuous  real-valued  functions  f and  g as 


a b a b 


f and  g over  the  interval  as  xs  b can  be  estimated  by  dividing 
the  interval  as  b into  a sufficiently  small  number  of  equal  length 
segments.  The  criteria  for  sufficiently  small  is  determined  by  the 
precision  of  the  information  desired.  For  further  information  on 
bandwidth  and  sampling  criteria,  see  "Introduction  to  Sonar  Technology", 
NAVSHIFS  0967-0C9-3010. 

Now  let  us  divide  f and  g over  as  xs  b into  n equal  segments.  Each 
segment  of  equal  length  will  have  a width  of  Ax,  Ax  - — — 


r 

re- 


consider the  following  sum: 
n-j 

E f(x1)g(x  )Ax 
i=0  1 1 J 


(2) 


Let  x = a and  x = b and  x.  be  the  i-th  samnle  of  f and  x.  .be 
o n i - i+j 

the  i+j-th  sample  of  g.  Let  Ax  be  a value  such  that  t = j Ax.  Then 

if  we  take  the  limit  as  Ax-0  of  expression  (2),  we  obtain: 

n-j  b-T 

lim  £ f(x  )g(x  )Ax  = j f(x)g(x+T)dx 

Ax-0  i=0  1 1+<3  a (3) 

which  is  the  numerator  of  equation  (l).  Thus,  expression  (2)  is  an 
estimate  of  the  numerator  of  equation  (l).  We  can  write  a similar 
equation  for  the  denominator.  Thus  we  have,  from  (l)  and  (3), 
b-T 

/ f(x)g(x+T)dx 

o(t)  = a 

b-T 

J |f(x)g(x+T)|dx 


*<*  . 

n-J 

E f(xi)g(x,  + ,)Ax 
- 1=0  J 

n-j 

to  E f(xi)g(x1+J 
= 1=0  J 
n-j 

to£  lf(xi)g(xi+i)l 

i=0 
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n-J 

I f(x.)g(x  ) 

o(j)  - 1-0  1 1 J (4) 

n-J 

I |f(x  )g(x  )1 
i«0  J 

This  last  expression  is  the  equation  used  in  the  running  program 
to  compute  the  cross-correlation  function,  e(j). 

B.  Clipped  Reference  Cross-Correlation 

The  CLIF'IK?  system  routine  is  designed  to  operate  in  conjunction 
with  formatted  digitized  data  from  the  FORMAT^-'  computer  program. 

F0RMA3$J  is  a computer  program  which  processed  and  formatted  the  input 
data  used  in  CLIPREF.  The  fonaat  of  this  data  is  described  in  Section 
VI,  C of  this  memorandum.  The  origin  of  this  data  is  sonar  acoustic 
data  which  was  recorded  at  sea  on  analog  taws.  A signal  of  pseudorandom 
noise  or  frequency  modulated  linear  sweep  was  sent  into  the  water.  As 
it  was  being  sent  into  the  water,  it  was  recorded  by  feedback  circuitry 
onto  analog  recording  tapes.  This  type  of  information  is  called  a 
reference  signal  or  simply  a reference.  After  the  signal  had  been 
transmitted  into  the  water j medium,  the  external  receiving  gates  were 
opened,  and  sonar  information  was  received  by  the  system  from  the  ocean 
itself.  This  information  was  also  recorded  on  analog  tapes  for  further 
analysis  and  is  called  an  echo  signal  or  more  simply,  a signal.  Later, 
at  th,.  laboratory,  these  echo  and  reference  blocks  of  data  containing 
analog  acoustic  data  were  digitized  by  the  FORMAT^-'  program  and  stored 
on  digital  magnetic  tape,  so  that  this  information  could  be  used  as 
data  input  to  analysis  programs  such  as  CIIPREF. 

Each  element  of  echo  or  reference  data  on  mag  tape  consists  of  1C 
bits  of  signed  information.  Inspection  of  Equation  ( 1 ) indicates 


that,  the  computation  of  a cross-correlation  value  involves  a great  deal 


of  multiplication.  To  multiply  two  signed  twelve  bit  numbers  requires 
from  35.2  to  112  microseconds  on  the  Univac  1206  computer.  Thus  to 
compute  just  one  element  of  this  sum  would  require  at  least  110  - 1^0 
microseconds.  Straight  computation  of  the  cross-correlation  function 
by  Equation  (4)  is  excessively  time  consuming  where  timing  is  important. 

However,  this  cross -cor relation  function  can  be  very  closely 
approximated  by  clipping  either  the  reference,  the  echo  or  both.  The 
effect  of  the  clipping  would  be  a greatly  simplified  computation  and 
therefore,  usually,  a greatly  reduced  computer  computation  time. 

Clipping  is  defined  to  be  that  process  by  which  a signal  or  element  of 
information  is  so  reduced  in  precision  as  to  be  limited  to  having  only 
two  possible  values,  relative  to  given  thresholds. 

If  we  arbitrarily  designate  these  two  values  as  zero  and  one  and 
the  reference  value  as  p,  then  we  shall  assign  a data  element  a value 
of  1 if  greater  than  or  equal  to  p and  a value  of  0 if  less  than  p. 

In  this  particular  case,  i.e.,  in  the  program  CLIPREF,  p=0  because  the 
information  desired  is  the  sign  of  reference  elements.  However,  the 
clipped  value  assigned  is  0 if  the  element  is  positive  and  1 if  the  value 
is  negative. 

In  the  CLIPREF  program,  we  are  clipping  the  reference.  Thus,  if  f 
is  the  echo  and  g is  the  reference,  Equation  (4)  reduces  to: 


c(j)  = i=0 


n-j 

E f(xt)g( 


W 


n-j 

E lf(xi)| 

i=0 


(5) 


where  g(*i+^)  = ±1.  With  Equation  (5)..  the  computation  of  a cross- 
correlation  value  can  be  done  in  the  computer  by  a comparison  of  signs 
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sequence  in  the  numerator  and  summation  and  a straight  summation  in 
the  denominator  rather  than  multiplication  in  both  numerator  and  denomi- 
nator. Equation  (5)  is  faster  than  Equation  (4);  however.,  the  technique 
actually  used  by  CLIPREF  is  much  faster  than  even  Equation  (5).  This 
is  disclosed  in  the  section  explaining  the  GENCORR  subroutine. 

C . Programming  Options 

There  are  basically  two  programming  options  available:  continuous 
mode  and  segmented  mode. 

1.  Continuous  Mode:  Continuous  mode  cross-correlation  consists 
of  the  cross-correlation  of  a short  signal  with  a long  signal  by  sliding 
the  shorter  signal  along  the  longer  signal  from  one  end  of  it  to  the  other. 
There  are  two  ways  in  which  this  can  happen  in  CLIPREF:  (l)  The  echo 
is  longer  than  the  reference.,  (2)  the  reference  is  longer  than  the  echo. 

We  wish  to  define  the  following  terms:  REFNUM  is  number  of 
reference  input  samples;  ECHNUM  is  the  number  of  echo  input  samples; 

JJ  is  the  number  of  output  cross-correlation  values.  Accordingly,  if 
REFNUM  > ECKNUM,  then  JJ  = REFNUM  - ECHNUM  + 1,  and  if  ECHNUM  > REFNUM, 
then  JJ  = ECHNUM  - REFNUM  + 1.  Let  E^  be  the  set  of  echo  signal 
samples  where  E^  is  the  i-th  sample  of  the  echo  signal  and  let  R^ 
be  the  set  of  reference  samples  where  R^  is  the  i-th  reference  sample. 

The  formulas  for  the  cross-correlation  functions  are  as  follows: 

Case  jtl:  REFNUM  > ECHNUM 
ECHNUM-1 

E R^.-E, 

o(j)  = 1=0  J x , J - 0,  1,  2 JJ  - 1 (6) 


ECHNUM-1 


. ■ , . 


I 


Case  4 2:  ECHNUM>REFNIIM 
REFNUM-1 

T.  <5  = 

i=0 

REFNUM-1 

t "W 

1=0 

2.  Segmented  Mode:  Segmented  cross-correlation  is  the  cross- 
correlation of  two  signals  of  equal  length  by  equi-length  segments. 
Each  segment  is  slid  a maximum  of  MM  points  in  each  direction  and  is 
correlated  against  its  counterpart.  QQ  is  the  number  of  segments 
into  which  the  input  data  shall  be  divided.  MM  is  the  maximum 
number  of  samples  that  the  echo  signal  is  slid  against  the 
stationary  reference.  The  number  of  output  samples  per  segment.  JJ_, 
is  equal  to  2MM+1.  The  total  number  of  output  samples  would  be 
QQ(2MM+l).  In  segmented  cross-correlation,  the  ECHO  segment  less 
MM  samples  at  each  end  is  moved  across  the  reference  from  the  left 
end  at  j = -MM  to  the  right  end  at  j = +MM.  Consider  the  following 
diagram  in  which  QQ  = 2. 


Ref: 


SEG 


Ech: 


MM 


VAL 


MM 


(8) 


In  the  CLIPREF  routine,  the  reference  and  echo  files  must  be 
equal  and  an  integral  multiple  of  QQ.  In  practice.,  if  they  are 
unequal,  the  longer  file  of  data  is  truncated  at  the  end  so  as  to 
equal  the  shorter  file.  Then  both  are  truncated  at  the  ends,  if 
necessary.,  to  make  them  an  integral  multiple  of  QQ.  SEG  is  the  number 
of  elements  per  segment  and  is  computed  by  SEG  = REFNUM/QQ.  VAL+1  is 
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is  the  number  of  elements  in  the  echo  segment  which  are  correlated 
with  the  reference.  VAL  = SEG  -(2MM+1).  Each  segment  is  computed 
by  the  following  equation: 

VAL+MM 

a(j)  = 1=MM  1+J  1 (9) 

VAL+MM 

V.  lEj 

i=MM 

where  i is  indexed  from  the  beginning  of  the  segment.  Inspection 
of  (9)  reveals  that  the  denominator  of  the  right  hand  side  is  a 
constant  with  respect  to  indexing  j,  and  hence  need  only  be  computed 
once. 

3.  Output  Options;  There  are  two  parts  to  the  output.  The 
first  part  consists  simply  of  a listing  in  order  of  computation  of 
all  the  valu's  computed.  The  values  are  simply  outputted  and  listed 
ten  to  a line  in  a table  similar  to  a log  table.  The  second  part 
of  the  output  consists  of  a graph  of  the- cross -correlation  function. 
The  graph  is  precise  to  0.02  or  to  2$.  However,  values  are  printed 
one  to  a line;  hence  the  second  ovtput  is  ten  times  longer  than  the 
first.  To  save  time,  the  second  output  can  be  ignored  if  so  desired 
by  setting  key  1.  If  the  second  output  is  desired,  tlrn  clear  key  1. 

Core  Option  for  Reference  Input;  Input  data  Is  read  in  to 
the  computer  memory  in  two  sets.  The  reference  data  is  read  first, 
then  the  echo  data.  Each  successive  run  requires  a different  set 
of  echo  data.  However,,  successive  runs  my  require  the  same  set 
of  reference  data.  To  save  time  on  data  input,  there  is  a core 
option  for  the  reference.  If  core  option  is  not  used,  a new  run 
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will  read  in  a new  reference  file.  If  core  option  is  selected., 
the  old  reference  stays  in  core  and  only  a new  echo  is  read  into 
the  computer.  Selection  of  core  option  is  discussed  below  under 
parameter  cards  preparation. 

II.  PROGRAM  USAGE 

A.  Card  Deck  Preparation 

1.  Job  Card:  Every  job  at  the  waterfront  computer 
laboratory.,  NUWC,  begins  with  a job  card.  The  job  card  should 
include  the  following  information:  (l)  the  code  making  use  of  the 
program.,  (2)  the  problem  number  funding  the  use  of  the  program,  (3) 
the  estimated  time  of  run.,  (4)  name  of  the  programmer  requesting 
the  run.,  and  (5)  the  date  on  which  the  run  was  submitted.  Job 
card  format  is  discussed  in  the  WFA  COMPUTER  CENTER  USERS  MANUAL. 

2.  LIB  Cards;  Two  library  routines  are  used  in  conjunction 
with  CLIPREF.,  PRINT  and  LILFLOP.  PRINT  is  a general  purpose  output 
routine  for  such  devices  as  high  speed  printer,  flexowriter  and 
paper  tape  punch.  LILFLOP  is  a general  purpose  routine  for 
performing  floating  point  operations.  It  is  used  by  CLIPREF 

to  format  the  print-out  of  decimal  numbers  for  the  PRINT  package. 

PRINT  is  based  at  14-0  and  runs  almost  to  2050.  LILFLOP  is  based 
at  2050  and  runs  almost  up  to  7000.  Library  cards  format  is  explained 
in  the  WFA  User's  Manual.  Load  PRINT  at  l40  and  LIT, FT, OP  at  2050. 

3.  L0D4q  and  TRA  Cards:  CLIPREF  is  based  at  7000  and  runs 
up  to  about  20500.  The  remainder  of  core  to  67777  is  used  as  a 
general  data  storage  area.  Currently  CLIPREF  is  on  File  7 of  tape 
WFA-268.  When  using  this  system.,  be  sure  to  mount  the  program 

tape  on  unit  M2.  After  loading  CLIPREF  and  the  two  library  routines, 
initiate  use  of  the  program  by  transferring  to  7000  with  a transfer 


card.  Format  of  L0D40  and  TRA  cards  is  also  described  in  WFA 
User's  Manual. 


B.  Parameter  Cards  Preparation 

The  parameter  cards  follow  the  transfer  card.  Each  distinct 
correlation  run  requires  two  parameter  cards.  There  is  no  limit 
to  the  number  of  runs  which  may  be  performed  after  the  transfer 
card.  A correlation  run  is  defined  to  be  the  cross-correlation  of 
two  files  of  data  as  specified  by  the  input  parameters  computed  in 
the  manner  specified  by  these  parameters.  A series  of  correlation 
runs  is  terminated  by  an  end  card. 

1.  First  Parameter  Card;  The  first  parameter  card  of  a 
correlation  run  contains  seven  parameters.  They  are  positive  decimal 
numbers  whose  values  do  not  exceed  five  digits.  They  are  named: 

(1)  REFFILE,  (f)  REFUNIT,  (3)  S IGF ILF,  ( h ) SIGUNIT,  (?)  QQ,  (6)  KM, 
and  (7)  HOPFLAG  and  their  positions  on  the  parameter  card  are  fields 
2 through  8 respectively,  each  parameter  occupying  one  field.  The 
most  significant  digit  of  each  parameter  must  be  punched  in  the 
left  most  digit  of  its  field;  otherwise  the  card  input  s/r  will 
consider  this  an  input  error.  WFA  User's  Manual  gives  a full 
description  of  card  fields  and  their  positions. 

REFFIIiE  is  the  file  on  mag  tape  on  which  the  reference  data  is 
located;  REFUNIT  contains  the  unit  number  of  this  tape.  SIGFILK 
contains  the  file  number  of  the  signal  echo  data  on  mag  tape  while 
SIGUNIT  has  the  unit  number  of  this  tape.  QQ  is  the  number  of 
segments  into  which  the  data  is  to  be  divided.  It  only  has 
significance  in  the  segmented  mode.  In  continuous  mode,  this  number 
is  zero.  MM  is  the  value  of  the  maximum  relative  time  delay.  Like 
QQ,  it  has  meaning  only  in  the  segmented  mode.  HOPFLAG  is  the 
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switch  indicator.  For  continuous  cross-correlation,  set  HOPFLAG  = Oj 
for  segmented  cross-correlation,  set  HOPFLAG  = 1. 

2.  Second  Parameter  Card:  The  second  parameter  card  of 
each  correlation  run  contains  the  time  code  limits  of  the  reference 
and  echo  files  of  the  data  to  be  selected.  The  four  parameters  on 

this  card  are:  (l)  initial  reference  time  code,  (2)  final  reference 

time  code,  (3)  initial  echo  time  code,  and  (4)  final  echo  time  code. 

Initial  reference  begins  in  field  2;  final  reference  begins  in  field 

5;  initial  echo  time  code  begins  in  field  8 and  the  final  echo  tine 
code  begins  in  field  11.  Each  time  code  consists  of  twelve  characters. 
There  are  two  hours  digits,  two  minutes  digits,  two  seconds  digits, 
and  three  milliseconds  digits.  Hours,  minutes,  seconds  and  milli- 
seconds fields  ai'e  sepai’ated  by  periods.  The  following  time  code 
example,  18.45. 08.100,  signifies  6:45  pm,  8-l/l0  seconds.  An 
example  of  these  parameter  cards  is  given  below. 

r T5745X87 1'OO TST45:-  f 3. 1 007  1 6 r^.'DSTJ  00 TBr4*n3.T0tr~ 

It 7 r«r»~rri[rrj  •<  nr  y lof.-  .ni .« 7V;  ~ v i>n> %r.  »r » * vjinW:.»s?~ ? r Vo]  >: « 7 ± 


oToT  

1 1 1 rrrr.’s » ■ i;  1 p>j. 

1 1 


13 3 13 "3" 


ID'  — TOO’I — 


72 
1 > 
) 1 

( 4 

1 » 

a 

Si 

I ? 

7 7 

I I 

1 1 

9 '9, 


■i  'I  lix  u n . ii :i 


ill  11  M 1'  i:  ~i  '*1.  • '*  imi'ii  . i j . ~ t u s?  i-n  l 1 m 1 -.  1 1 p m ino] . 11' - j • 


IS* 

mi 


2 2 2 2| 
* > > 

3 133 


4 4 4 4 

i «! 


tin 


2 2 2 2 
ii)i 

3 3 3 3 


4 4 4 4 

Sill 


• 10  n 

11*1 


"O'Bjfl  0 0 Cju'o  0 0{0  o o'c}j  c o c]c  C 0 Oil)  0 0 cTl 


0 0 0 OjOOOOjO'O  0 G,0  C 0 C|0  C 0 OjO  0 0 CjG  0 o c;o  0 C 00  G oT', , , 

1 ti  m ii  ii'ii  ii  i)  k .i  mm.iiii  n »,»»  n t:]i>  i<  u kh  » ;i  sjjn  .1  uujixi  11  si'ojosi  sijsisisicln  tin 

I*  1 t il  1 


I I 11 


tilt  11  1 1 11  1 11 


It  >1111  11111111111 


2 2 2 2!2  2 2 212  2 2 2 2 2 2 2’2  2 2 21?  2 2 2 2 2 2 2: 2 2 2 7 2 2 2 7 2 2 2 212  2 2 2 2 2 2 2l2  7 2 2 2 2 2 2 7 


1111-1111 


• 10  II  l.'i-l  11  IS  IS 


1 3 3 


4 4 4 

1 10  n >:] 


5 5 5 515  5 5 5-5  55  5 


6 6 6 6| 
s ; s • 
7 17  7 


3 3 3 


4 4 4 4 

is  r •!  ii 


I I I I I 


1 - --  --  - | - 1 — — — 1 

11  is  ii  10:11  min  :s  n 11  in  is  n ::  n u h i:  1 mu;  ;>  II  Ii.-y  u si  1 11  n n v->i  si  <s  ;.'n  11  n«,i  P 1 1 ts 

3 3 3 3 3 ' 3 3 3 | 3 3;3  3 3 3;3  3 3 1*3  3 3 3.3  3 3 1:3  3 3 ill  3 3 3,3  3 3 3‘3  3 3 3|3  3 3 3 3 


4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4'4  4 4 4 1 4 4 4 4 4 4 4 4 4 4 4[4  4 4 4't  4 4 4|4  4 4 4, 1 . 

»>  il  >1  ;4«  a *J  .*<  ?•  r|.'i  M :i  J.’  » ?» ?» ’ ■« ;ii;  n i:  i;  u n « it  • • *»  e y «•  : ss  mu’*-  si n *•  «;  tj r ,* 

5 5 5 5’5  5 5 5 5 5 5 5 0 5 5 5-555  j's  5 5 5 5 5 5 55  5 5 5!5  5 5 5 5 5 5 5*5  5 5 5f5  5 5 5.j  5 5 5 5 : 

| | I j I r QfNtiAi.  • :o  n*.i8i  III 

6 6 4 6 6 8 6 6,8  S 8 6 8 8 6 6 5 E S E G £ 6 E|6  6 6 £;S  fS"?  7i  8 $ 6 * 9 8 3 8 t't  f!  8 8 5 8 * St  S S 3 6 6 3 f • 

j*  7 vt  .3  ii  iii'j  n «j  13*17  it  -■»  :c*.m  » *<  r u y r :i  nr  , - ' «n  «t  • 'n??v  i.!>m  u.-Vjif*  . 

m Hii  nil  n ill  i nil  inli  i n.ii  n\n  inn  hi  i n n n\ii  i il;  j 1 ill  iiilin  n . 

il!  I 

I I 8 S S 8 3 8 S 3 8 o’l  1 i 3|S  S S ! 9 ! I S S I S 3, S S 1 8 S 1 S o 3 3 I ! S S 3 S I J I 3 3 I S i S 1 1 l-i  3 S I S i 

3 » > » 1 30  « mi  :yj«  :?  ji  !» « jj  ji  j»  ii  j: )« *14- » »» ••  i «•«  om  ji  »*' i uh>'  pouii 

9 9 9 9:9  3 9 S 3 9 3 ?!$  9 3 S 9 S 9 «;9  9 9 5 9 3 9 9 8 3 3 9'3  3 3 p 9 9 3 3 3 9 15  3 5 5 5|3  3 3 3 S 9 3 3 2 « 9 5 J 5 S 9 < 

I t 1 «i)  I 7 ill  10  II1MI  M tmllMI  II  joljl  Jt /Is.’IIJ  II  JSI1M  JlJ8'r  )U»«.  «•  «MMI  4' •!  it  • 


118  8 

1 J > I 


\>\iV  V:'|4J»W  ».  •»*„ 


11 


5.  Core  Option  for  Reference;  To  select  the  CORE  option 


for  reference  as  described  above,  delete  tooth  tiraecodes  for  the 
reference  and  punch  the  word  CORE  in  Field  2. 

C.  END  Cards 

An  END  card  terminates  a sequence  of  correlation  runs.  It 
returns  the  program  to  MONITOR  control. 

III.  PROGRAM  OUTLINE:  GENERAL  DESCRIPTION 

The  executive  routine  CLIPCORR  calls  upon  the  following  major 
subroutines  to  perform  the  task  of  the  program:  (l)  CARDINPUT, 

(2)  TRANTIME,  (3)  READPACK,  (U)  RDUNPACK,  (5)  COMPCORR,  (6)  PRINTOUT. 

A.  CARDINPUT 

CARDINPUT  reads  in  the  two  parameter  cards,  translates  the 
information  and  stores  the  values. 

B.  TRANTIME 

TRANTIME  translates  the  echo  and  reference  initial  and  final 
time  code  tables  to  milliseconds.  Inputs  are  the  nine  word  tables 
REFLIMKl,  REFLIMHI,  ECHLIMLO  and  ECKLIMHI.  Outputs  are  the  cells 
labeled  MILIREFLO,  MILIREFHI,  MILIECHLO,  and  MILIECHHI. 

C.  READPACK 

The  subroutine  READPACK  reads  in  reference  data  from  magnetic 
tape  as  specified  by  the  time  codes,  clips  it  retaining  only  the 
sign  information  of  each  element  and  packs  it  into  core  memory  30 
elements  to  the  word.  Inputs  to  READPACK  are  mag  tape,  FILE  and 
UNIT  information,  and  the  reference  time  limits.  Output  is  into  the 
general  core  area  and  the  table  limits  word,  REFTAB. 

D.  RDUNPACK 

This  procedure  reads  and  unpacks  an  echo  signal  from  mag  tape 
into  the  core  area.  Inputs  to  RDUNPACK  are  FIIiR  and  UNIT  information, 


mag  tape  and  the  echo  time  code  limits.  The  data  are  read  in, 
unpacked  and  sign  extended  from  twelve  bits  to  fifteen  bits  and 
stored  in  the  output  table  AREA  as  half  words.  The  word  AREA  contains 
the  upper  and  lower  limits  of  the  AREA  table. 

E.  COMPCORR 

This  routine  is  the  heart'  of  the  CLIPREF  system.  It  does 
all  the  housekeeping,  updating  and  so  on  for  the  various  correlation 
computations.  It  initializes  and  executes  the  subroutine,  GENCORR, 
which  performs  cross -correlations  as  specified  by  inputs.  Inputs 
to  COMPCORR  are  such  quantities  as  QQ,  MM,  HOPFLAG,  FKFLAG,  ECHNUM, 
and  REFNUM.  The  output  of  COMPCORR  is  a number  of  cross-correlation 
data  written  on  a scratch  file  of  the  program  tape. 

F.  PRINTOUT 

This  routine  prints  the  output  of  a correlation  run. 

Input  is  from  the  program  mag  tape.  Output  is  on  the  high  speed 
printer.  There  are  two  outputs:  (l)  an  output  of  data  and  index 
values  and,  (2)  a graphical  output  of  the  data.  The  segmented 
mode  of  the  graphical  output  prints  out  as  many  graphs  as  segments 
with  the  time  codes  of  each  segment  on  the  header  of  the  graph. 

IV.  CARDINPUT 

A.  Translation  of  First  Card  Parameters 

The  parameter  card  is  first  read  into  core  using  the  CARD  s/r 
in  MONITOR.  The  parameter  cells  REFFILE,  RFFUNIT , ECHFIIE,  ECHUNIT, 
QQ,  MM,  and  HOPFLAG  are  in  order.  They  are  cleared.  The  information 
in  CARPRF,  the  input  buffer  used  by  CARD,  is  translated  word  by  word 
from  XS-3  codes  to  octal  numbers  and  stored.  Each  parameter  is 
assumed  to  be  an  unsigned  decimal  integer.  Each  parametric  field 
is  unpacked  character  by  character  into  a nine  word  table  called  PACK. 
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The  unpacking  loop  then  tests  the  characters  and  stops  when  a blank 
or  an  illegal  character  is  reached.  An  error  condition  is  generated 

, I 

when  either  of  the  following  conditions  exist:  (l)  The  first 
character  of  a field  is  blank  or  (2)  an  illegal  character  code  is  in 
the  field,  i.e. , a non-numeric  character  other  than  a blank.  If 
such  an  ei*ror  is  generated,  the  following  message  is  typed  on  the 
flex:  FIRST  PARAMETER  CARD  ILLEGAL.  The  computer  than  U- stops  at 
the  call  to  CARPINFUT  in  the  EXEC.  All  the  operator  need  do  then  is 
punch  a correct  card,  stick  it  at  the  bottom  of  the  hopper,  and 
hit  high  speed.  If  all  is  well,  the  parameters  are  properly  stored 
and  the  second  parameter  card  sequence  is  initiated. 

B.  Translation  of  Second  Card  Parameters 

Immediately  after  the  second  parameter  card  is  read  into 
core,  the  third  word  of  CARP  Kb’  is  checked  to  see  if  it  contains  the 
word  CORE.  If  it  does,  the  subroutine  Jumps  to  the  translation  of 
the  echo  time  codes,  skipping  translation  of  the  reference  time 
codes.  Each  time  code  parameter  is  unpacked  by  the  same  loop.  The 
loop  first  unpacks  each  parameter  into  PACK,  then  it  transfers  PACK 
to  the  appropriate  9 word  table.  These  output  tables  are  REFLIMLO, 

REFLIMHI,  ECHLIMLO,  and  ECHLIMHI.  As  each  parameter  digit  is  unpacked, 
it  is  checked  for  numeric  value.  If  any  of  the  digits  checked  are 
non-numeric,  the  following  message  is  printed  on  the  flexowriter; 

SECOND  PARAMETER  CARD  ILLEGAL.  The  program  then  stops  at  the  call 
to  CARPINPUr  in  the  EXEC.  The  card  may  then  be  corrected  by  the 
operator;  however,  both  cards  must  be  read  to  restart  the  Job.  If 
all  is  well,  after  the  parameter  tables  are  filled,  CARP INPUT  transfers 
control  to  the  BiXEC. 
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C.  Recognition  of  End  Card;  Exit  Condition 


i 


I 


If  the  first  card  read  in  by  CARDINPUT  is  an  end  card., 
control  is  transferred  back  to  the  MONITOR. 

V.  TRANIIME 

A.  Specific  Explanation 

TRANTIME  translates  echo  and  reference  initial  and  final 
time  tables  to  milliseconds.  A subroutine  named  PACKTIME  is  used 
to  do  this.  The  four  input  tables  each  containing  nine  time  code 
digits.,  REFLIMLO,  REFLIMHI , ECHLIMLO,  and  ECHLIMHI  are  translated 
to  milliseconds  and  stored  in  cells  MILIREFLO,  MILIREFHI.,  MILIECHLO ., 
and  MILIECHHI  respectively. 

PACKTIME  translates  timeburst  coded  time  to  millisecond 
time.  Input  is  the  nine  word  TIMER  table;  output  is  the  cell 
MILITIME . 

B.  Compute  REFNUM  and  ECKNUM 

TRANTIME  computes  REFNUM  and  ECHNUM  by  the  following  formulae: 

REFNUM  = (MILIREFHI  - MILIREFLO ) (25 )/2 

ECHNUM  = (MILIECHHI  - MILIECHLO ) (2 5 )/2 
TRANTIME  then  returns  control  to  the  EXEC. 

VI.  READPACK 

A.  Specific  Information 

This  procedure  reads  in  reference  data  and  packs  it  into 
core.  The  contents  of  REFFILE  and  REFUNIT  are  stored  into  FILE 
and  UNIT  and  A is  cleared  in  order  to  initialize  POSTIME.  POSTIME 
positions  the  specified  tape  unit  to  the  specified  cell,  searches 
for  the  first  time  code  of  the  data  and  moves  it  to  TIMER.  READPACK 
then  executes  PACKTIME  and  stores  MILITIME  into  REFSTART.  REFSTART 
is  compared  with  MILIREFLO . If  MILIREFLO  is  smaller  than  REFSTART, 


the  followirg  error  is  printed  onto  the  flexowriter:  INITIAL  T3MEC0DE 
NOT  IN  FILE.  If  all  is  well,  the  subroutine  SHOVEMIN  is  initialized 
and  executed.  SHOVEMIN  packs  the  reference  data  into  core.  After' 
execution  of  it,  the  reference  data  table  limits  cell  REFTAB,  is 
filled  and  READFACK  returns  control  to  the  exec. 

B.  POSTIME 

POSTIME  has  two  functions.  On  the  first  call  of  a data 
input  sequence.,  it  positions  a tape,  spaces  past  the  header  record, 
reads  the  first  data  record  and  transfers  the  first  time  burst  to 
TIMER.  On  subsequent  calls , it  merely  reads  in  an  additional  record 
of  data.  Inputs  are  the  cells  FILE  and  UNIT.  Outputs  are  the 
tables  TIMER  and  DIG. 

The  A-register  is  the  switch  for  POSTIME.  If  A is  zero, 
POSTIME  assumes  a first  call.  If  A is  non-zero,  it  merely  reads 
another  record.  If  on  an  initial  call  to  POSTIME  the  formatted 
digitized  data  is  searched  and  no  time  code  is  found,  the  following 
message  is  printed  on  the  flexowriter:  INITIAL  TIMECODE  MISSING. 

The  job  is  then  aborted  with  o core  dump  being  printed  on  the  high 
speed  printer. 

Input  data  from  mag  tape  is  read  into  table  PIG.  PIG  is  a 
1200  decimal  word  table  designed  to  use  and  handle  input  data  only. 

If,  when  reading  a record  of  data,  TAPE  reads  an  end  of  file,  FOSTIME 
executes  an  abnormal  exit.  If  POSTIME  is  executed  by  READPACK,  the 
following  message  is  typed  on  the  flexowriter:  FIRST  REFERENCE  RECORD 
MISSING. 

C.  Mag  Tape  Input  Pata  Format 

A complete  section  on  quantity  or  input  data  is  called  a 
file.  Each  file  of  data  consists  of  a set  of  records . There  is  no 
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set  number  of  records  to  a file;  however,  each  record  Is  exactly 
1200  words  long.  The  first  record  of  each  file  is  a header  record 
and  contains  general  Information  about  the  data  to  follow.  Subsequent 
records  contain  the  formatted  digitized  data  itself.  In  the  last 
record  of  a file,  the  end  of  the  data  is  marked  by  the  filling  of 
the  rest  of  the  table  with  negative  zeroes.  Following  this  record 
is  an  end  of  file  marker. 

The  data  word  format  is  as  follows: 
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IT 


CHAN  A 


CHAN  B 


TIMECODE 


12  bits  12  bits  6 bits 

The  upper  24  bits  of  the  word  are  divided  into  two  fields  of 
12  bits  each  called  CHAN  A and  CHAN  B.  The  lower  six  bits  contain 
the  time  code  information.  There  are  two  modes  in  which  data  can 
be  packed  into  these  records  and  stored  on  tape,  dual  mode  ard 
single  mode.  In  dual  mode,  two  separate  records  of  data  are  stored., 
one  in  CHAN  A and  the  other  in  CHAN  B.  In  single  mode  storage,  the 
information  is  stored  alternately  in  CHAN  A and  CHAN  B,  two  samples 
to  the  word  with  there  being  only  one  channel  of  data.  All  input 
data  for  CLIPREF  is  single  mode  data. 

Every  hundred  samples,  a 13  word  time  code  burst  is  generated. 

In  single  mode  format,  this  burst  occurs  every  50  words.  The  time 
code  burst  contains  13  binary  coded  decimal  numbers.  These  numbers 
are  in  bits  0-3  in  the  TIMECODE  field.  Bit  4 is  the  time  burst 
marker,  and  when  lit,  indicates  the  first  word  of  a time  burst. 

The  first  three  codes  of  a burst  are  milliseconds;  the  next  two 
are  seconds;  the  two  following  are  minutes  digits;  the  next  two  are 


IT 
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hours.  The  next  three  are  the  run  number;  and  the  last  digit  is  a 
code  number.  A more  complete  description  of  the  data  format  is 
described  in  the  TRG  report,  TRG  0G3-TN-66  entitled:  Sea  Trial 
Data  Acquisition  and  Analysis  Equipment. 

Although  a time  code  burst  occurs  every  50  samples  throughout 
the  input  data,  the  input  routines  of  CLIPREF  translate  only  the 
first  time  code  burst  of  each  file  of  data  and  ignore  all  the  rest. 

D.  Core  Data  Format 

1.  GENERAL  Table  Design:  Library  routines  run  from  address 
lUO  to  7000.  The  CLIPREF  program  itself  starts  at  7000.  The 
MONGOOSE  MONITOR  runs  from  70000  on  up.  CLIPREF  with  all  its 
instructions  and  fixed  length  tables  runs  from  7000  to  about  20500. 
This  leaves  a general  data  storage  area  from  approximately  20500 

to  70000.  This  area  is  called  the  GENERAL  table,  and  it  is  here  that 
reference  and  echo  data  are  stored.  Reference  and  echo  data  are 
not  fixed  length  data,  hence,  they  are  stored  in  core  in  relative 
table  areas. 

Reference  data  is  packed  in  first,  the  starting  address 
being  the  starting  address  of  GENERAL.  When  packing  is  completed, 
the  upper  and  lower  address  limits  of  the  table  are  stored  in  the 
cell  labeled  REFTAB,  the  upper  limit  in  the  upper  half  and  the  lower 
limit  being  stored  in  the  lower  half.  The  echo  data  is  packed  in 
core  immediately  above  the  reference  data.  After  storage  into  core, 
its  table  limits  are  stored  into  AREA  in  the  same  manner. 

2.  REFTAB  Table  Design:  The  reference  table  is  packed 
containing  only  sign  information.  The  data  are  packed  30  bits  to 
the  word,  starting  from  the  left  going  to  the  right.  A set  bit 
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indicates  a negative  sign;  a cleared  bit  indicates  a positive  sign. 

In  order  to  perform  any  correlation  at  all,  all  the  reference  data 
must  be  read  into  core.  Since  the  available  cells  number  20000. 
this  limits  the  reference  data  input  to  about  600.000. 

E.  SHOVEMIN 

This  procedure  packs  one  bit  sign  data  from  mag  tape  30 
bits  to  the  word  and  stores  it  as  specified  by  Bl,  which  contains 
the  address  to  store  the  packed  value.  Inputs  are  Bl,  BB5  and  NN. 

BB5  contains  the  initial  number  of  values  to  skip  before  starting 
the  packing  cf  data.  NN  contains  the  number  of  data  points  to 
pick  up.  Outputs  are  REFTAB  and  the  table  to  which  it  refers. 

The  subroutine  first  spaces  ahead  the  proper  number  of 
values,  making  use  of  P0STIME  to  do  this.  If  P0STIME  reads  an  end 
of  file,  an  error  is  generated  and  the  following  message  is  typed 
on  the  flexowriter:  INITIAL  TIMEC0DE  NOT  IN  FILE.  After  spacing 
forward  to  the  correct  starting  point,  SHOVEMIN  starts  packing  the 
data  and  stores  them  in  the  reference  table.  As  each  data  word 
is  pulled  out  of  table  DIG,  it  is  tested  for  a value  of  negative 
zero.  If  a negative  zero  is  found.,  the  following  error  is  typed 
on  the  flexowriter:  LIMErs  EXCEED  DATA  ON  FILE.  If  the  length 
of  the  reference  data  exceeds  the  core  limits  (about  600,000),  an 
error  is  generated  and  this  message  is  typed  on  the  flexowriter: 

CORE  OVERFLOW.  If  all  is  well,  after  packing  in  the  data.,  SHOVEMIN 

. 

returns  control  to  CLIPCORR  and  the  system  proceeds  to  read  the 

\ 

echo  data. 

VII.  RDUNPACK 


This  procedure  reads  and  unpacks  an  echo  signal  from  mag  tape 
into  the  relative  table  AREA.  Its  operation  is  comparable  to  RFAPPACK. 


ECHFILE  and  ECHUNIT  are  stored  into  FILE  and  UNIT.  The  time  code 


parameters  are  converted  to  milliseconds.  PASSFORD  is  used  to 
space  to  the  beginning  of  the  signal  from  the  beginning  of  the 
file.  The  echo  date  is  packed  in  a different  format,  but  other  than 
that,  packing  and  storage  is  similar  to  READFACK. 

A.  PASSFORD 

This  procedure  passes  forward  in  order  to  reach  the  first 
cell  of  data.  BB5  is  used  as  the  number  of  samples  to  skip. 

POSTIME  is  used  to  actually  pass  records  forward. 

B.  AREA  Table  Design 

The  AREA  table  starts  directly  after  the  reference  table 
and  ends  at  address  67777.  Each  word  of  AREA  contains  two  samples 
of  data  at  most.  The  input  data  are  signed  12-bit  numbers.  When 
packed  into  AREA,  they  are  sign-extended  to  half  words.  Data  is 
packed  into  AREA  sequentially  in  the  upper  half  of  the  word  f u 
the  lower  limit  to  as  high  as  67777.  When  the  upper  half  words 
have  been  filled,  then  data  packing  starts  in  the  lower  half  words 
from  the  lower  limit  to  67777.  In  the  course  of  such  packing, 
one  of  three  conditions  will  arise:  (l)  data  input  will  not  completely 
fill  the  upper  half  words  of  the  table,  (2)  data  input  will  fill  the 
upper  half  but  not  the  lower  half  words  of  the  table,  (3)  data  will 
completely  fill  the  table  and  not  all  the  echo  can  be  read  into  core. 

If  (l)  occurs,  then  the  upper  limit  of  AREA  is  set  where  it  stopped 
and  program  is  exited.  If  (2)  occurs,  then  the  upper  limit  of 
AREA  is  set  to  67777  and  subroutine  is  exited.  If  (3)  occurs,  then 
the  upper  limit  of  AREA  is  set  to  67777.  a flag  is  set  carrying  the 
information  that  not  all  the  data  could  be  read  in.  and  RDUNFACK  is 
exited. 
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C.  Output  Function  Words 


The  following  variables  are  also  computed  by  RDUNPACK  in 
the  course  of  reading  the  echo  signal  data  into  core.  They  are 
subsequently  used  by  COMPCORR  to  determine  appropriate  courses  of 
action. 

1.  FKFLAG:  FKFLAG  is  the  packed  echo  signal  flag.  If 
FKFLAG  = 0,  then  all  echo  data  selected  has  been  packed  into  core. 
FKFLAG  = 1 indicates  that  not  all  signal  data  was  packed  into  core. 

2.  DIGPT:  DIGPT  is  the  last  input  word  referenced.  It 
contains  the  index  number  of  table  DIG  of  the  last  signal  datum 
stored  into  AREA.  If,  upon  subsequent  execution  of  COMPCORR, 

it  is  desired  to  read  an  additional  amount  of  echo  signal  data, 

DIGPT  indicates  the  point  from  which  the  last  sample  had  been  read. 

5-  DCTR:  DCTR  is  the  echo  data  counter.  It  contains  the 
number  of  echo  samples  which  have  been  read  in  a run  regardless  of 
the  number  actually  present  in  core.  DCTR  is  usually  equal  to  NAR, 
but  if  AREA  has  been  updated  to  input  more  data,  it  will  contain 
the  number  of  total  data. 

4.  ICTR:  ICTR  is  the  initial  data  item  counter.  It 
contains  the  index  number  of  the  first  echo  sample  currently  in  core. 
It  is  cleared  by  RDUNPACK,  and  remains  zero  until  AREA  is  updated 
and  subsequent  data  is  read,  e.g.,  suppose  ICTR  = 0 and  AREA  is 
updated  by  reading  J00  more  samples.  The  whole  table  must  be  slid 
down  300  samples  with  the  bottom  300  deleted  from  the  table.  ICTR 

is  then  incremented  by  300. 

5.  NAR:  NAR  is  the  number  of  AREA  items.  It  contains  the 
total  number  of  echo  signal  samples  presently  in  core.  The  following 
relationship  holds  among  these  last  three  quantities: 
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NAR  + ICTK  = DCTR 


r 


Nr 
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VIII.  COMPCORR 

The  COMPCORR  routine  with  i+s  attendant  subroutines  constitutes 
the  heart  of  the  CLIPREF  system.  The  other  parts  of  the  system 
perform  merely  input,  output  and  housekeeping  functions.  COMPCORR 
does  the  actual  cross-correlation  computations. 

Upon  entry  to  COMPCORR,  the  program  tape  is  positioned  to  the 
scratch  file  and  the  cross-correlation  output  table.  CCOUT,  is  cleared. 
HOFFLAG  is  interrogated  and  control  is  given  over  to  the  appropriate 
mode,  either  conti nuous  or  segmented. 

A.  Continuous  Mode 

The  definition  and  selection  of  continuous  mode  cross - 
correlation  has  already  been  discussed.  Four  distinct  conditions 
can  arise  in  continuous  correlation:  (l)  All  signal  data  are  in 
core  and  the  signal  is  longer  than  the  reference,  (2)  all  signal 
data  are  in  core  and  the  reference  is  longer  than  the  signal,  (?)  not 
all  signal  data  are  in  core  and  the  signal  is  longer  than  the  reference, 
and  (4)  not  all  signal  data  are  in  core  and  the  reference  is  longer 
than  the  signal. 

1.  ECHNlftbREFNUM,  PKFLAO  ==  0:  This  is  probably  the  case 
which  is  used  most  often.  All  correlation  is  computed  with  one  call 
to  GENCORR  with  no  updating.  The  correlation  output  values  are 
computed  by  Eq(7). 

2.  ECHNUM~>RKFNliM , rKFLAG  ~ 3:  The  correlation  output 
values  are  still  computed  by  Eq('f)  in  this  case,  but  now  the  AREA 
table  must  be  updated  from  time  to  time  since  the  signal  data  are 
not  all  in. 
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3.  REFNlflfoECHNUM.  FKFLAG  = 0:  All  correlation  values  are 


computed  in  one  call  to  GENCORR  with  no  updating.  Eq(6)  is  used 
to  make  the  computations. 

k.  RKFNUM>EC HNUM , FKFIAfl  = 1;  This  condition  is  impossible 
because,  no  correlation  can  be  computed  by  Eq(6)  unless  all  echo 
data  are  in  core.  An  error  is  generated  and  the  following  message 
is  typed  out:  CORRELATION  IMPOSSIBLE  ECHO  FILE  TOO  LONG. 

B.  Segmented  Mode 

Similarly,  the  definition  and  selection  of  the  segmented 
mode  has  aLready  been  discussed.  However,  in  this  mode,  the 
reference  and  the  signal  must  be  made  to  be  the  same  length.  Each 
segment  to  be  correlated  must  satisfy  the  following  condition: 

SEGVMM  + 1.  If  this  is  not  the  case,  then  an  error  is  generated 
and  the  following  message  is  typed:  MM  TOO  LARGE  FOR  SEGMENT  LENGTH. 
Three  cases  can  arise  in  segmented  cross-correlation:  (l)  all 
data  in  core,  (2)  some  data  not  in  core , but  at  least  one  whole 
segment  is  in  core,  and  (3)  some  data  not  in  core  and  part  of  first 
segment  not  in  core. 

Case  (l):  Segmented  cross-correlation  performed  without 

updati ng. 

Case  (2 ) : Segmented  cross-correlation  performed  with  updating. 

Case  (3):  The  segment  is  too  long  to  perform  a cross- 
correlation. An  error  is  generated  and  the  following  message  is 
typed:  CORRELATION  IMPOSSIBLE,  SEGMENT  TOO  LONG. 

C.  GTFLAG  Modes 

GTFLAG  stands  for  greater  than  flag.  It  has  the  value  of  0 
if  RKFNUM>ECHNUM  and  the  value  of  1 otherwise.  GTFLAG  is  an  input 
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of  GENCORR  and  essentially  tells  it  whether  to  J- index  through 
the  echo  data  or  through  the  reference  data.  In  continuous  mode 
GTFLAG  is  1 for  cases  1 and  2 and  0 for  cases  3 and  4.  GTFLAG 
is  always  0 for  segmented  mode  cases. 

D.  FKFLAG  Modes 

HCFLAG  indicates  whether  or  not  all  the  data  have  been  read. 

If  all  echo  data  are  in,  correlations  may  always  be  made.  If  not, 
there  are  cases  in  which  correlations  are  impossible.  In  continuous 
mode,  this  occurs  when  the  reference  is  longer  than  the  echo.  In 
segmented  mode,  this  occurs  when  the  segment  length  is  longer  than 
the  signal  in  core. 

E.  UPDATE  S/R 

This  procedure  updates  table  AREA  by  an  amount  equal  to  ELO 
and  reads  in  more  data  as  available.  ELO  is  mainly  an  input  for 
GENCORR  and  is  the  initial  echo  signal  index.  UPDATE  is  called  only 
after  GENCORR  has  been  executed.  When  UPDATE  is  called  EI0  + VAT. 
has  exceeded  NAR,  the  number  of  echo  values  in  AREA.  One  of  three 
conditions  may  occur  when  UFDATE  is  called:  (l)  ELO>NN,  (2)  ELOs 
NN,  and  (3)  abort  condition.  After  GENCORR  has  been  executed,  NN 
contains  the  number  of  whole  words  in  the  table  AREA. 

1.  Case  1;  ELO>NN:  In  this  case  ELO  samples  on  the  bottom 
will  be  deleted.  Samples  in  the  lower  half  of  the  table,  starting 
with  the  ELOth  sample  will  be  moved  to  the  beginning.  When  all 
samples  in  the  table  have  been  moved,  the  rest  of  the  table  is  cleared 
and  the  input  of  additional  data  is  initiated.  Data  input  at  this 
stage  is  similar  to  the  data  input  processing  of  RDUNPACK. 


2.  Case  2;  ELO<flJN:  ELO  number  of  samples  at  the  bottom 


of  the  table  are  deleted.  Samples  in  the  upper  half  of  the  table 
move  up  ELO  positions.  At  the  bottom  of  the  upper  half,,  are 
inserted  samples  from  the  top  of  the  lower  half.  Then  the  lower 
half  table  is  updated.  New  samples  are  packed  in  at  the  bottom 
of  the  lower  half. 

In  both  cases  1 and  2,  ICTR  will  be  incremented  by 
ELO.  If  all, data  is  read  in,  FKFLAG  will  be  cleared  and  NAR  will 
be  updated  so  as  to  correspond  to  the  actual  number  of  samples 
currently  in  AREA. 

3.  Abort  Condition:  An  abort  condition  will  occur  in  the 
data  input  stage  of  UPDATE  under  the  same  conditions  that  would  cause 
an  error  generation  in  the  data  input  of  RDUNPACK.  The  same  action 
will  be  taken  as  in  RDUNPACK. 

F.  GENCORR 

COMPCORR  essentially  does  the  housekeeping  and  initialisation 
for  the  cross-correlation.  GENCORR.,  the  largest  subroutine  in 
CLIPRKF computes  the  correlation  values  and  writes  them  on  mag 

tape. 

1.  Indexing  and  Input  Parameters:  Inputs  are  RLO,  ELO., 

VAL,  NAR,  JJ,  and  GTFLAG,  Indices  used  are  Bl,  B2,  B3,  BB4,  and  B?. 
RLO  is  the  initial  reference  index.  It  contains  the  index  number 
of  the  first  datum  to  be  pulled  out  of  the  reference  table.  ELO 
is  the  initial  echo  index  and  contains  the  index  of  first  datum 
to  be  pulled  out  of  the  echo  signal  table.  VAL  is  the  segment  of  the 
correlated  interval.  It  contains  the  number  of  correlation  pairs 
that  are  to  be  compared.  JJ  is  the  correlation  index  shift.  It 
contains  the  number  of  values  and  shifts  to  be  computed  in  any  call  to 


GENCORR. 
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B1  is  the  reference  input  index.  It  contains  the  address 
of  the  word  in  REFTAB  containing  data  to  be  correlated.  B2  is  the 
echo  input  index.  It  contains  the  address  of  a word  in  AREA  to  be 
correlated  with  the  reference.  B3  is  the  value  counter  index.  It 
counts  the  number  of  pairs  being  correlated  and  is  indexed 
through  VAL.  BBU  is  the  shift  counter  index.  It  counts  the  number 
of  shifts  of  the  data  and  is  indexed  through  JJ.  B5  is  the  output 
counter  index.  It  counts  the  number  of  data  computed,  and  when  the 
output  table,  CCOUT,  is  filled,  CCOUT  is  written  on  tape,  it  is 
cleared,  and  B5  is  reinitiated. 

In  the  initial  housekeeping,  NN  is  computed.  NN  is  the 
number  of  samples  in  the  upper  half  of  AREA.  It  is  a value  which 
is  used  constantly  in  the  indexing  of  the  correlation.  Next  the 
denominator  of  the  right  sides  of  Eq(6 ) or  Eq(7)  is  computed. 
Immediately  after  this  computation  of  the  absolute  sum,  which  is  done 
only  once  per  call  to  GENCORR,  the  general  indexing  cycles  begin. 

The  general  indexing  cycle,  which  must  be  done  for  each 
correlation  shift,  is  rather  complex.  Reference  data  which  may 
come  from  any  bit  of  a 30  bit  word  must  be  matched  against  signal 
data  which  may  come  from  an  upper  word  or  from  a lower  half  word. 
Reference  data  may  start  anywhere  in  a 30  data  input  word  and  may 
end  anywhere  in  a 30  data  input  word.  A signal  data  segment  may 
lie  entirely  in  the  upper  half  of  AREA,  it  may  lie  entirely  in  the 
lower  half  of  the  table,  or  it  may  spill  over,  some  data  in  the 
upper  half  and  some  in  the  lower  half. 
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The  actual  correlation  Is  done  by  two  min  cross-correlation 
sequences.,  each  of  which  will  correlate  JO  elements  at  a time. 

One  sequence  correlates  reference  data  with  upper  half  signal 
data;  the  other  correlates  it  with  lower  half  signal  data.  The 
index  matching  is  divided  into  six  phases;  three  for  the  upper 
sequence  and  three  for  the  lower  sequence.  The  variables  for  these 
phases  are  01,  02,  0J,  0i+,  05  and  06. 

2.  Upper  Half  Table;  The  general  indexing  cycle  must  first 
determine  whether  signal  values  will  be  pulled  from  upper  half  only, 
lower  half  only,  or  from  both.  If  the  signal  must  be  pulled  from 
upper  half  only  or  from  both  halves,  the  upper  half  sequence  must 
be  initiated.  The  variables  01,  02.,  and  0J  must  be  computed  for 
this  sequence. 

01  is  the  index  number  of  the  first  reference  value  in  a 
reference  data  word  corresponding  to  RLO  in  the  reference  data. 

To  compute  01,  divide  RLO  by  ^0  and  put  the  remainder  in  01 . The 
quotient  is  added  to  the  lower  half  of  REP TAR  and  stored  in  Rl. 

Phase  1 picks  up  the  initial  fraction  of  reference  data  in  the  upper 
half  sequence.  Phase  2 correlates  whole  reference  words,  phase  J 
picks  up  the  final  fractional  reference  word  for  correlation.  02 
contains  the  number  of  times  phase  2 is  repeated.  Add  01  to  the 
number  of  upper  half  values  and  subtract  JO;  then  divide  by  JO. 

The  quotient  is  02;  the  remainder  is  0%  0J  contains  the  number  < 

of  reference  values  to  be  picked  in  phase  J.  After  computation  of 

these  parameters,  the  indexing  cycle  computes  the  lower  sequence 

parameters  if  signal  data  carry  over  to  the  lower  half;  otherwise 

control  is  transferred  to  phase  1 correction.  A special  case  occurs 


when  the  reference  data  are  less  than  30  and  all  in  the  same  word. 

9(1  is  computed,  but  control  is  transferred  directly  to  phase  3, 
bypassing  phases  1 and  2. 

3.  Indexing  for  Lower  Half  Table:  Indexing  for  lower 
half  correlation  sequence,  which  involves  the  computation  of  04, 

05.  and  06  is  greatly  similar  to  the  indexing  of  the  upper  half 
sequence.  04,  05,  and  06  functionally  correspond  to  01,  02  and  03 
respectively. 

4.  Updating  Sums;  Consider  Eq(6)  and  Eq(7).  If  correlation 
is  computed  by  Eq(6),  then  the  denominator  need  not  be  updated  at 

the  end  of  a correlation.  At  the  end  of  each  correlation,  GENCORR 
checks  GTFLAG.  If  GTFLAG  = 0,  we  are  using  Eq(6)  and  we  increment 
RLO  and  jump  to  the  general  index  cycle.  Otherwise  we  are  using 
Eq(7)  and  must  update  the  denominator.  GENCORR  simply  subtracts 
the  absolute  value  of  the  first  value  and  adds  the  absolute  value 
of  the  next  new  value,  according  to  the  following  relation: 


s,.,  ■ Si  - |!W  + "WvM,1  “here  Si 


EL0-1+VAL 

V IE  I 
i = ELO  1 i' 


5.  Method  of  Computation:  Eq(6)  and  Eq  (7)  ai'e  actually 
too  slow  to  be  used  in  the  computer  for  the  actual  correlation 
computation  and  have  been  given  here  for  heuristic  purposes.  Consider 
the  numerator  of  Eq  (7): 


£ RE.,  J = 0,  1,  ....  JJ-1 
i=l  1 J 


Each  element  of  the  sum  must  be  either  plus  or  minus  |E^+j|.  We 
define  the  following  functions: 
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r E., . when  R.E.^,  is  positive 

Ei+J(+)  = i+J  1 i+J 

' 0 otherwise 


p -E...  when  is  negative 

Ei+i(.)  * i+J  1 i+J 

' 0 otherwise 


Accordingly  (l)  becomes: 


i=!  RiEi+J  = Jx  V)(+)  “ V>(-) 


We  also  have  the  following  relation  for  the  denominator  of  Eq  (7): 


• VAL 

S4  = v |E 


j = i=!  ,Ei+J!  = ^ Ei+J(+)  + Ei+j(-) 


iSi  Ei+J(+)  = SJ  ‘ ^ Ei+J(-) 


By  substitution,  Eq  (7)  becomes: 


o(j)  - Jl  ~ = SJ  ~ 1?!  E^J(-)  (1j) 


A similar  relation  holds  for  Eq  (6): 


a(J)= 


S - 2 y E 

r>  L 


h i(-) 

i=l  v ‘ , J = 0,  1 JJ-1 

S 


Where 


0 if  R.  ,E.  positive 
F (-)  = i+J  i 

iv  ' -E^  otherwise 


and  where 


VAIj 

S°  ‘ ili  |Ei' 


***** 
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Eq(l3)  and  Eq(lU)  are  the  equations  actually  used  in  the 
computation  of  the  correlation  values.  The  matching  of  data 
pairs  thus  involves  only  the  computation  of: 

VAL  VAL 

£ Vk-)  or  & V) 

6.  Output  Format:  The  computed  correlation  values  are 
stored  in  a whole  word  in  CCOUT  scaled  up  28  hits. 

DC.  PRINTOUT 

PRINTOUT  is  the  output  subroutine.  It  reads  in  from  mag  tape 
and  prints  it  out  as  previously  explained.  The  output  is  in  two 
parts,  the  correlation  value  print-out  and  the  cross-correlation 
graph  print -out. 

A.  Correlation  Value  Print-Out 

For  this  output  the  program  generates  a page  feed,  prints 
the  title  "clipped  reference  cross-correlation" , then  the  date 
below  it,  a table  header,  and  lastly,  the  correlation  values  to 
five  decimal  places  in  order  as  they  appear  on  tape. 

B.  Cross-Correlation  Graph  Print-Out 

For  this  output  PRINTOUT  generates  a page  feed,  prints  the 
title  header  and  the  date  below.  All  the  input  parameters  follow. 
Included  is  a three  line  explanation  of  the  graph.  Lastly  follows 
the  graph,  itself. 

1.  Key  1 Option;  This  output  can  be  deleted  by  setting  Key  1. 

2.  Continuous  Mode  Graph:  The  continuous  mode  graph's 
index  starts  at  zero. 

3*  Segmented  Mode  Graphs;  Segmented  mode  graph  may  have 
more  than  one  segment.  In  each  segment,  the  index  starts  at  -MM. 


The  time  interval  for  echo  and  reference  is  divided  into  QQ  equal 
intervals.  The  time  interval  for  each  such  segment  is  printed  at 
the  top  of  each  segment's  graph. 

C.  UNPACKTIME 

UNPACKTIME  is  a subroutine  used  by  PRINTOUT  to  set  up  the 
output  of  the  time  codes.  It  translates  millisecond  time  to 
timeburst  coded  time.  Input  is  MILITIME  and  output  is  the  TIMER 
table . 

D.  SETEMUP 

SETEMUP  is  a subroutine  used  by  PRINTOUT  to  set  up  the  output 
of  the  time  codes.  It  translates  the  time  codes  as  stored  in  the 
TIMER  table  and  generates  a string  of  literals  for  output  on  the 
printer  in  conjunction  with  the  PRINT  package.  Input  is  TIMER; 
output  is  STRING2. 

E.  GRAFH 

This  procedure  is  used  by  PRINTOUT  to  print  a cross-correlation 
graph.  Inputs  are  CCOUT  and  FONUM.  PONTJM  is  the  number  of  lines 
to  output.  GRAPH  sets  up  horizontal  lines  every  tenth  column  on 
every  line.  Every  tenth  line,  vertical  bars  are  set  in  the  printer 
output  table,  GRAF.  GRAFH  uses  the  subroutines,.  SENDWO  to  send  a 
line. 

1.  Output  Format:  101  columns  starting  at  column  19  make 
up  a line.  Correlation  values  are  indicated  by  an  asterisk.  All 
output  values  are  rounded  off  to  the  nearest  0.02. 

2.  SUBGRAF:  This  procedure  stores  a print  out  value  for 
output  in  table  GRAF. 


3.  SENDWO:  This  procedure  sends  the  printer  19  blanks, 
then  it  unpacks  the  codes  in  GRAF  and  sends  them  to  the  printer. 
After  this,  it  prints  a line  without  a line  feed. 


F.  PRTDATE 

This  procedure  is  used  by  PRINTOUT  to  print  the  current 

date. 

G.  SETTIME 

This  procedure  computes  a time  code  and  sets  up  a string 
of  literals  to  print  it  out.  Inputs  are  REFLIMLO  and  Bl.  Output 
is  STRING2;  subroutines  used  are  UNPACKTIME  and  SETEMUP. 

X.  TIMING 

A.  Data  Input 

Translation  and  read-in  of  parameter  cards  takes  less  than 
two  seconds.  The  reading  and  packing  time  of  input  data  is  usually 
a matter  taking  less  than  20  to  30  seconds.  If  the  data  is  a long 
ways  down  the  tape,  the  tape  may  take  awhile  to  wind  down  to  it,  but 
usually  input  data  time  considerations  are  trivial. 

B.  Cross-Correlation  Time  Estimation 

The  time  for  cross-correlation  depends  upon  four  main 
factors:  (l)  the  average  indexing  time,  (2)  the  average  time  for 

matching  one  pair  of  samples,  (3)  the  length  of  the  correlation 
segment,  and  (4)  the  number  of  correlation  output  values . 

The  general  index  cycle  time  varies  from  a low  value  of 
839-2  microseconds  to  1904.8  microseconds,  depending  upon  the  path 
taken  and  GTFLAG  mode.  The  average  index  time  cycle  for  GTFLAG  = 0 
is  1267.9  microseconds  whereas  indexing  the  j-shift  through  the  echo 
signal  takes  an  average  of  1475-9  microseconds,  about  208  microseconds 
longer.  Let  I be  the  average  indexing  time,  which  we  will  take  as 
1475.9  microseconds. 
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Let  T be  the  average  time  for  matching  one  pair  of  samples, 

S be  the  number  of  correlation  pairs,  and  J be  the  number  of 
correlation  output  values.  The  time  to  compute  one  output  value 
would  be  ST+I,  hence  the  total  time,  t,  would  be  t = j(St+I). 

The  numerical  value  for  T is  t = 52.4  microseconds.  For  S 
small,  the  indexing  time  I is  a large  percentage  of  the  time  required 
for  correlation.  As  S gets  larger,  the  effect  of  I on  the 
correlation  time  gets  smaller  and  smaller.  Consider  the  following 
table: 


Sample  Size 

Time  (usee) 

$ Indexing 

10 

2000 

73-7# 

30 

3048 

48.4$ 

100 

6716 

22.0$ 

300 

17196 

8.57$ 

1000 

53876 

2.74$ 

3000 

158,676 

0.93$ 

10000 

525,476 

0.28$ 

30000 

1,573,476 

0.09$ 

C.  CC  Value  Print-Out 

Time 

To  print  one  line 

of  the  ten  value  per  line 

output  takes 

about  0.250  sec/line  +2$. 

D.  Cross-Correlation  Graph  Print-Out  Time 

To  print  one  line 

of  the  graph  takes  about 

0.0945 

sec/line  +2$. 

SUMMARY 

The  CLIPREF  system  was 

based  on  the  experience  gathered  from 

the  production  of  the  CLIPCORR  programming  system. 

In  the  CLEPREF 

system,  it  was  important  to  speed  up  the  computations;  This  was 
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done  at  the  expense  of  considerable  core.  Computation  elements 
which  would  have  to  be  repeated  millions  of  tines  in  a correlation 
run  were  strung  out  thirty  times  per  loop  instead  of  being  repeated 
merely  once  per  loop.  The  result  was  highly  successful  in  terms 
of  correlation  computation  speed.  The  value  of  T was  reduced 
from  about  80  microseconds  to  52. ^ microseconds. 

Additional  information  is  included  in  the  appendices. 

Appendix  A is  a list  of  possible  error  printouts  and  their 
significance.  Appendix  B is  the  flow  charts  of  the  main 
subroutines  of  the  system.  Appendix  C contains  examples  of  the 
outputs.  Both  outputs  to  two  continuous  correlation  runs  are 
shown.  The  first  graph  shows  the  data  points  unconnected.  In 
the  second  graph  we  have  connected  the  data  points  to  show  the 
sinusoidal  nature  of  the  data  which  we  were  correlating. 
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APPENDIX  A 


Error  Print-Outs 

1.  FIRST  PARAMETER  CARD  ILLEGAL  - Illegal  format  on  first 
parameter  card.  Program  will  4-STOP  at  the  beginning  of  the 
card  input  sequence.  Correct  the  parameter  card  and  place  it  at 
the  beginning  of  the  deck.  Restart  the  program  by  hitting  high 
speed. 

2.  SECOND  PARAMETER  CARD  ILLEGAL  - Illegal  format  on  second 
parameter  card.  Program  will  4-STOP  at  the  beginning  of  the 
card  input  sequence.  Correct  the  parameter  card  and  place  it 
after  the  first  card  at  the  beginning  of  the  deck. 

3.  INITIAL  TIMECODE  MISSING  - Indicates  that  no  timecode  has 
been  found  on  the  first  record  of  data  read  into  core.  This 
usually  indicates  a bad  record  of  data  on  mag  tape.  Frogram  will 
give  a core  dump  and  terminate  the  job. 

4.  CORE  OVERFLOW  - Indicates  that  the  amount  of  reference  data 
has  exceeded  the  general  area  core  limits  assigned  to  it.  Program 
will  give  a core  dump  and  terminate  the  job. 

5.  LIMITS  EXCEED  DATA  ON  FILE  - Indicates  that  one  of  the  parametri 
upper  timecode  limits  of  the  input  data  is  greater  than  the  last 
timecode  of  the  corresponding  data  file  on  mag  tape.  Therefore, 
that  data  cannot  be  read  and  unpacked  into  core.  The  program  will 
Jump  to  0 and  Fault. 

6.  FIRST  REFERENCE  RECORD  MISSING  - Indicates  that  the  lower 
reference  timecode  parameter  limit  cannot  be  found  on  the  reference 
input  tape.  The  program  will  jump  to  0 and  Eau.lt. 
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7-  INITIAL  TIMECODE  NOT  IN  FILE  - Indicates  that  the  initial 
reference  timecode  is  less  than  the  first  timecode  on  the 
reference  data  tape,  and  hence  the  timecode  is  not  on  the  file. 
Program  will  Fault  and  jump  to  0. 

8.  FIRST  SIGNAL  RECORD  MISSING  - Indicates  that  the  lower  signal 
timecode  exceeds  the  last  timecode  on  the  signal  input  tape. 

Similar  to  FIRST  REFERENCE  RECORD  MISSING.  Program  will  jump  to 
0 and  Fault. 

9.  CORRELATION  EXCEED  DATA  AVAILABLE  - Indicates  that  an  error 
has  been  made  in  a call  to  GENCORR.  The  number  of  correlations 
called  for  has  exceeded  the  amount  of  signal  data  in  core 
necessary  to  make  the  computations. 

10.  CORRELATION  IMPOSSIBLE,  ECHO  FILE  TOO  LONG  - This  error 
occurs  when  the  reference  is  longer  than  the  echo  data  and  all  the 
echo  data  cannot  be  read  into  core  because  of  lack  of  core 
storage.  Thus  the  correlation  cannot  be  performed  because  the 
echo  file  is  too  long.  The  program  jumps  to  0 and  Faults. 

11.  MM  TOO  LONG  FOR  SEGMENT  LENGTH  - This  error  print-out  can 
occur  only  when  the  segmented  correlation  mode  has  been  selected. 

It  indicates  that  MM,  the  number  of  correlation  shifts  is  larger 
than  half  the  segment  length.  Correlation  therefore  cannot  be 
performed.  Program  jumps  to  0 and  Faults. 

12.  CORRELATION  IMPOSSIBLE , SEGMENT  TOO  LONG  - This  error  can 
occur  only  when  the  segmented  correlation  mode  has  been  selected. 

It  indicates  that  the  signal  segment  is  too  long  for  all  of  it  to 
be  read  into  coi’e.  The  program  ignores  the  run  causing  this  error. 


jumps  to  the  beginning  of  the  executive  routine  and  begins  executing 


the  next  correlation  run  request. 
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APPENDIX  C 


CLIPCORR  S/R 


CLEAR  GENERAL 
CCEE  STORAGE 


CAP.Dirpir 


f I, 


COMPCORR 


ENTER 


/ POSITION  \ 
' PROCRAM  TAPE  > 
TO  SCRATCH  FILE 


CLEAR 


CCP.RELATIOM 
OUTPUT  TABU 


CONTINUOUS 

MODE 


ECHNUI:  > 


->GTF1AG 


( GTFLAO  = 

0 

\ YES 

- 

J 

1 

1 


COMPCORR 
Pg  2 


COMPCORR 


CORRELATE 

ifposs  im: 

ECHO  FILE 
TOC  LORO 


INITIAL! 

GEKCCEE 


FAULT 


' CORRELATE  JJ  > 
VALUES  INDEXING 


LFNCTI 

IR’LTIF 


COMPCORR 

Pg  5 


INITIAI.IZE 

GEKCCRR 


IS  UPDATE 
NECESSARY 


/READ  ID  TCP!' 
ECHO  SAIPIJ'S 
AS  AVAILAPLE 


\ UPDATE 


*■>•  V .D  j i.  I o 

caunrED 


PRINTOUT 


ENTER 


f POSITION  \ 
OUTPUT  TAPE 
THE  BEGINNING 
OF  THE  DATA 


READ  ONE  > 
RECORD  OF 
OUTPUT  DATA 


PRINT  CROSS' 
CCEEELATION 
HEADER 


PRINTOUT 

Pg  3 


PRINT 

DATE 


PP.TDATE 


PRINT  SIGNAD 
AND  REFERENCE 
FILE  NUMBERS 


PRINT 


T PRINT  SIGNAL' 
AMD  REFERENCE 
TAPE  IT’ IT 
NT  TERRS 


PRINT  QQ  ' 
AND  MM  VALUE; 


CONTINUOUS 


PRINTOUT 
Pg  6 


' SET  UP  FOR 
PRINT  CUT  OF 
FINAL  ECHO 
T I NEC  ODE 


PRINT  ECHO 
TIMECODFS 


PRINT 


COMPUTE  LONER 
REFERENCE  TCI 


UITACKTi; 


SET  IT  FOP 
PRINT  CUT  CF 
INITIAL  REFER] 
THiECODE 


CONFUTE  UPPER 
REFERENCE  TINECOi 


TRA?!SL\TE  > 
:.M,  EFFERENT' 
E FN  SEC  ('NR, 

* J u'ArUlJi'.tJ 


